February 4, 2008 20:56 Proceedings Trim Size: 9in x 6in dfest0116 



in 
o 

O ■ EINSTEIN- iETHER THEORY* 

(N 

CHRISTOPHER ELING t , TED JACOBSON 1 ''*, AND DAVID MATTINGLY § 

t Dept. of Physics, University of Maryland, College Park, MD 20742-4111, USA 
* Institut d'Astrophysique de Paris, 98bis Bvd. Arago, 75014 Paris, France 

§ Dept of Physlcs ' Unwersity of CaUfomia > Dams ' Cahfomm ' 95616 ' USA 

> 



We review the status of "Einstein- jEther theory", a generally covariant theory of 
gravity coupled to a dynamical, unit timelike vector field that breaks local Lorentz 
symmetry. Aspects of waves, stars, black holes, and cosmology are discussed, to- 
gether with theoretical and observational constraints. Open questions are stressed. 
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Q-i 1. Introduction 

^Jy Could there be an aether after all and we have just not yet noticed it? By 

an "aether" of course we do not mean to suggest a mechanical medium 
whose deformations correspond to electromagnetic fields, but rather a lo- 
■ cally preferred state of rest at each point of spacetime, determined by some 

hitherto unknown physics. Such a frame would not be determined by a 
circumstance such as the moon's gravitational tidal field, or the thermal 
cosmic microwave background radiation, but rather would be inherent and 
unavoidable. Considerations of quantum gravity have in multiple ways led 
to this question, and it has also been asked in the context of cosmology, 
where various puzzles hint that perhaps something basic is missing in the 
standard relativistic framework. 

Lorentz symmetry violation by preferred frame effects has been much 
studied in non-gravitational physics, and is currently receiving attention 
as a possible window on quantum gravity. 1 But what about gravity itself? 
General relativity is based on local Lorentz invariance, so if the latter is 
violated what becomes of the former? It is hard to imagine, both philosoph- 
ically and technically, how we could possibly give up general covariance, the 
deep symmetry finally grasped through Einstein's long struggle. Thus the 



*Based on a talk given by T. Jacobson at the Deserfest. 
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question that interests us here is whether a generally covariant effective 
field theory with a preferred frame could describe nature. a 

The simplest description of such a frame would appear to be via a 
scalar field T, a cosmic time function, which has been proposed in various 
contexts. 3 ' 4 ' 5,6 The gradient T a , if timclike, defines a preferred rest frame, 
and one can envision dynamics that would force it to be everywhere time- 
like. But while a scalar field is simplest, the norm of the gradient \T a \ is 
"extra information" , which has nothing to do with specifying a frame per 
se but rather specifics the rate of a particular cosmic clock. It may be that 
Nature provides such a clock; we just wish to point out that the clock rate 
is extra information. Constraining the gradient to have fixed norm is not a 
viable option since, as explained in section 3, this would lead inevitably to 
caustics where T. a b diverges. Another noteworthy feature of using a scalar 
is that, by construction, the 4- velocity of the preferred frame is necessar- 
ily hypersurface-orthogonal, i.e. orthogonal to the surfaces of constant T. 
Again, perhaps this is the way Nature works, but it is a presumption not in- 
herent in the notion of a local preferred frame determined by microphysics. 

The alternative discussed in this paper is to describe the preferred frame 
by a vector field constrained kinematically to be timelike and of unit norm, 
which we call the (Ether field u a . Such a field is specified by three in- 
dependent parameters at each point, and generally couples via covariant 
derivatives, so the theory is far more complicated than that of a scalar time 
function. It is instinctive to worry about ghost modes given a vector field 
without gauge invariance. However the unit constraint on the vector ren- 
ders it an unfamiliar beast. All variations of the vector are spacelike, since 
they connect two points on the unit hyperboloid, so ghosts need not arise. 

There is a sparse history of studies of unit vector fields coupled to 
gravity. 7,8 ' 9 ' 10 ' 12 ' 13,14 Here we focus on the particular approach and results 
in which we have been involved. 15 ' 16 ' 17 ' 18 ' 19 ' 20 We begin with the action 
principle that defines the theory, and then discuss a Maxwell-like special 
case, linearized waves, PPN parameters, energy, stars and black holes, and 
cosmology. 

2. Einstein-aether action principle 

In the spirit of effective field theory, we consider a derivative expansion of 
the action for the metric g a b and aether u a . The most general action that 



a More general sorts of Lorentz violation in the gravitational sector are examined in Ref.2. 
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is diffcomorphism-invariant and quadratic in derivatives is 

S = Jd^x^^R + K ab mn V a u m V b u n + \{u a u a - 1)) (1) 

where 

K ab mn = c ig ab g mn + c 2 8 a m 8 b n + cz8 a J b m + c iU a u b g mn . (2) 

The coefficients 01,2,3,4 are dimensionless constants, R is the Ricci scalar, 
and A is a Lagrange multiplier that enforces the unit constraint. The metric 

signature is (H ), and units are chosen such that the speed of light 

defined by the metric g a b is unity. The constant G is related to the Newton 
constant Gn by a Q-dependent rescaling to be discussed below. Other 
than the signature choice we use the conventions of Rcf. 21. The possible 
term R a bU a u b is proportional to the difference of the c 2 and C3 terms via 
integration by parts, hence has been omitted. We have also omitted any 
matter coupling since we are interested here in the dynamics of the metric- 
aether sector in vacuum. Note that since the covariant derivative of u a 
involves the Levi-Civita connection, which involves first derivatives of the 
metric, the aether part of the action in effect contributes also to the metric 
kinetic terms. We call the theory with this action Einstein- (ether theory, 
and abbreviate using "iE-theory" . 

Another way to express the theory is using a tetrad e a A rather than 
the metric, where A is a Lorentz index. Then the aether can be specified 
as u a — u A e a Al with a unit Lorentz 4- vector u A satisfying the constraint 
Vabu a u b — 1, where ijab is the fixed Minkowski metric. This decouples the 
normalization condition on u A from the dynamical metric. The Lagrange 
density is then of the form K^qD a u A D\,vP , where D a is the Lorentz- 
covariant derivative involving the spin-connection u>^ D , and K AB is a linear 
combination of the four terms g ab r]AB, e A e b B , e B e b Al and u c u D e^e^riAB- 
This theory has a local Lorentz invariance, which can be used to set the 
components of u A to (1,0,0,0). That produces the form of the theory as 
presented by Gaspcrini. 7 One can also use a Palatini formalism, in which 
the spin connection is treated as an independent variable to be determined 
via its field equation. In this case the spin connection has torsion, because 
of the coupling to u A . If the solution ui^ D (e,u) is substituted back into 
the action, one returns to the tetrad form, but with different coefficients 
for each of the four terms in K^g. The relation between these coefficients 
and the original ones has not yet been worked out. 

Given a metric and a unit vector field, there is a one parameter family 
of metrics that can be constructed (aside from simple rescalings). When 
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expressed in terms of a different metric in this family, the action changes, 
but only insofar as the values of the Cj in (2) are concerned. More precisely, 
consider a field redefinition of the form 

g'ab = 9ab +{B- l)u a U b , u' a = Cu", (3) 

with C = [1 + (B - l)u 2 ] -1 / 2 (where u a := g am u m and u 2 = g mn u m u n ). 
Lorentz signature of both g ab and g' ab requires B > 0. The coefficient C is 
chosen such that g' ab u' a u' h = g a bu a u b , so the unit constraint is unchanged, 
hence in the action we can put simply C = B^ 1 / 2 . The action (1) for 
(g' ab , u' a ) takes the same form as a functional of (g a b, u a ), but with different 
values of the constants Cj. The general relation between the Cj and the c[ 
has recently been worked out by Foster 11 . His results reveal, for example, 
that one can arrange for c\ + c 3 = by choosing B = 1 — c[ — c' 3 (provided 
(c[ + Cg) < 1 for Lorentz signature). A special case previously worked out 
by Barbero and Villasehor 10 shows that the ae-theory is equivalent via field 
redefinition to GR when the parameters satisfy 13 c\ + c 4 = 0, ci + c 2 + c 3 = 0, 
and C3 = ±yjc\{ci — 2). Lorentz signature implies c\ < 0. Note that if one 
first makes a field redefinition such that c\ + C3 = 0, then the Barbero- 
Villasenor result reduces to the statement that the theory is equivalent to 
GR only if all coefficients vanish. Hence with the c\ + c 3 = description 
we ensure that non-zero coefficients always represent true deviations from 
GR. 

The field equation from varying the aether in the action (1) takes the 
form 

V a J a m - c 4 u a \7 m u a = Xu m (4) 

where 

J a m := K ab mn V b u n , (5) 

and u n = u b X7 b u n . The field equation from varying the metric in the action 
(1) together with a matter action takes the form 

G ab = T^ ] +^GTX ttc \ (6) 



b This corrects our earlier statement of the equivalence in Ref. 18. We thank B.Z. Foster 
for pointing out this error. 



February 4, 2008 20:56 Proceedings Trim Size: 9in x 6in 



dfcst0116 



5 

where the aether stress tensor is given by 18 

= V m (J (a m Uh) - J m [a u b) - J (ab) u m ) 

+C1 [(V m Ua)(V m U b ) - (V a U m )(V b U m )] 

+c 4 ii a u b 

+ K(V m J mn ) - c 4 u 2 ] u a u b 

-^9abL u . (7) 

The constraint has been used in (7) to eliminate the term that arises from 
varying \p^g in the constraint term in (1), and A has been eliminated by 
solving for it via the contraction of the aether field equation (4) with u a . 
The notation L u — —K ab mn VaU m VbU n is the aether lagrangian. 

Our goal is to determine the theoretical and observational constraints 
on the parameters Cj, and to identify phenomena whose observation could 
reveal the existence of the aether field. For such phenomena one can look 
at post-Newtonian effects, gravitational and aether waves, and cosmology. 



3. Maxwell-like simplified theory 

Before considering the general, rather complicated, theory it makes good 
sense to ask if there is a simplification that might serve at least as a decent 
starting point. A great simplification occurs with the choice ci + C3 = 
and C2 = c 4 = 0, so that the connections drop out of the aether terms in 
(1). The aether part of the Lagrange density then reduces to 

2 Cl u [aM u^+\(u 2 -l). (8) 

This theory was studied long ago by Nambu 22 in a flat space context, ft 
is almost equivalent to Einstein-Maxwell theory in a gauge with u 2 = 1. 
The difference is that one equation is missing, since the action need only be 
stationary under those variations of u that preserve u 2 = 1 . The missing 
equation is an initial value constraint equation, Gauss' law. If the current 
to which u a is coupled is conserved, then Gauss' law holds at all times if it 
holds at one time. 

This theory coupled to dynamical gravity was first examined in Ref. 8, 
and further studied extensively in Ref. 15 and Ref. 9. In Ref. 15 it was 
shown to be equivalent to Einstein-Maxwell theory coupled to a charged 
dust, restricted to the sector in which there exists a gauge such that the 
vector potential is proportional to the 4- velocity of the dust, i.e. the aether 
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field u a . The charge-to-mass ratio of the dust is given by c (8-itG/ci) 1 / 2 . 
The extremal value corresponds to c\ — 2. A number of results were estab- 
lished concerning static, spherically symmetric solutions, black holes, and 
linearized solutions. 

This case is appealing due to its simplicity, however a serious flaw was 
noticed: solutions can have "shocks" or caustics beyond which the evolution 
of the aether cannot be extended. In particular 15 , consider the aether config- 
urations that can be written as the gradient of a scalar u a — T a , so that the 
Maxwell-like "field strength" tensor U[ ab ] vanishes and u a is orthogonal to 
the surfaces of constant T. Then the field equations reduce to the vacuum 
Einstein equation, together with the vanishing of the Lagrange multiplier 
A and the statement that the gradient of T is a unit vector, T, a T' a = 1. 
Then u a is necessarily a geodesic: 

u a \7 a u b = u a \7 b u a = 0. (9) 

The first equality holds since u b is a gradient, and the second holds since it 
is a unit vector. If we launch geodesies orthogonal to an initial surface of 
constant T, they will generically cross after some finite proper time. Where 
they cross there is no well-defined value of u a , and the derivative W a u b is 
singular. These are the shock discontinuities. A different demonstration of 
the existence of shocks appears in Ref. 9. 

We note in passing that the preceding demonstration of shocks applies 
in a very different context, namely the version of fc-essence 4 recently called 
"ghost condensation" 5,6 . d This is the theory of a scalar field <j) with La- 
grangian density of the form P(X), where X = <l).a4>' a , where P(X) has 
a minimum at some value c. Among the solutions to the field equations 
is the special class which have X = c everywhere. The above argument 
shows that generically such configurations have caustics where <j)- a b is sin- 
gular. In the cosmological setting, Hubble friction drives all solutions to 
the minimum X = c. As far as we know it is an open question whether 
the diverging gradients of the X = c solutions is reflected in the generic 
cosmological solution. 

Another flaw with the Maxwell-like case is that it admits negative energy 
configurations, as shown by Clayton 9 using the Hamiltonian formalism in 
the decoupling limit where gravity is neglected. 6 The existence of negative 
energy configurations in this case is related to the fact that the Lagrange 

c The convention with Maxwell Lagrangian given by — F 2 /16n is adopted here. 

d For a discussion of caustics in a more general fc-essence scenario see Ref. 23. 

°The argument in Ref. 9 has a minor flaw, but the conclusion is correct. The negative 
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multiplier A can be negative, so in the charged dust interpretation the mass 
density is negative. 

Before returning to the general class of Lagrangians we note that one 
might also consider the theory where the restriction on the norm of u a 
is enforced not rigidly by a constraint but rather by a potential energy 
term V(u a u a ) in the action. This approach was discussed by Kostelecky 
and Samuel 8 , and more recently explored by Bjorken 24 , Moffat 25 , and 
Gripaios 26 . It has an additional, massive, mode, which should be checked 
for a possible wrong sign of the kinetic energy. 

4. Waves 

The spectrum of linearized waves is important for several reasons. First, 
it can be used to constrain the theory a priori, by rejecting values of the 
parameters Ci for which waves carry negative energy or for which there 
are exponentially growing modes. Second, wave phenomena can be used 
to place observational constraints on the parameters, using radiation from 
compact objects such as the binary pulsar, as well as cosmological pertur- 
bations. 

The spectrum of linearized waves around a flat spacetime background 
was worked out for the general theory defined by the action (1) in Ref. 19. f 
The wave modes in a de Sitter background were found in Ref. 14 (for 
c 4 = 0), which also studied the metric perturbations in inflation interacting 
with the vector as well as a scalar inflaton. 

Here we summarize the results for the modes around flat space. Since 
the aether has three degrees of freedom, the total number of coupled metric- 
aether modes is five. There are two purely gravitational (spin-2) modes, 
two transverse aether (spin-1) modes in which the aether vector wiggles 
perpendicular to the propagation direction, and one longitudinal or "trace" 
(spin-0) mode. The waves all have a frequency that is proportional to the 
wave vector. Hence they are "massless" and have fixed speeds. The speeds 
for the different types of modes are all different, and each mode has a 
particular polarization type. 



energy configuration described there is a time-independent pure gradient Uj = di<j>(x). 
This initial data with vanishing time derivative Uj^ = indeed has negative energy, 
however the equation of motion implies that the time derivative does not remain zero 
(unless Ui = 0). 

f The Maxwell-like special case was previously treated in Ref. 15, and the case with only 
ci non-zero was treated in Ref. 16. In Ref. 14 the modes were found in the small c; limit 
where the aether decouples from the metric (c/. section 6.2.2). 
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Table 1 gives the speeds and polarizations for the spin-2, spin-1, and 
spin-0 modes, in that order. The metric and aether have been expanded 
as g a b = Vab + h a b and u a = u a + v a , where i] a b is the Minkowski metric 
and u a is the constant background value that has components (1, 0, 0, 0) in 
the coordinate system adopted. The gauge conditions hoi = and Vij = 
are imposed, where i stands for the spatial components. The propagation 
direction corresponds to i = 3, and 7 = 1,2 labels the transverse directions. 
The notation C123 stands for ci + c 2 + c 3 , etc, and s is the wave speed. 

Table 1. Wave mode speeds and polarizations. 

squared speed polarization 

1/(1 -C13) h 12 , hu = -h 22 

(ci - \c\ + |c§)/ci 4 (l - C13) ft/3 = [ci 3 /s(l - c vi )]vi 

C123(2 - C14)/C14(l - Ci 3 )(2 + C13 + 3c 2 ) hoo = — 2l>0 

hn = h 22 = -ci4?Jo 
ft-33 = (2ci4/ci23)(l + C2)-»0 



4.1. Stability 

The squared speed refers to the squared ratio of frequency to wave- vector, 
so if it is negative for real wave- vectors the frequency is imaginary, implying 
the existence of exponentially growing modes. 6 The requirement that no 
such modes exist restricts the parameters of the theory. 11 For a small 
compared to unity this requirement reduces to the conditions C1/C14 > 
for the transverse vector-metric modes and C123/C14 > for the trace mode. 

Lim argued 14 that one should additionally demand that the modes prop- 
agate subluminally (relative to the metric g a b)- Although there is nothing 
wrong with local superluminal propagation in a Lorentz-violating theory, 
he pointed out that the vector field (in an inhomogeneous background) 
might tilt in such a way as to allow energy on such locally superluminal 
paths to flow around a closed spacetime curve. It is not clear to us that 
it is really necessary to impose this extra demand, since even in general 
relativity the classical field equations do not forbid the formation of closed 
timclike curves, around which relativistic fields could propagate. In any 

s The factor s in the polarization hj- A of the transverse aether mode implies that when 
s 2 < there is a 7r/2 phase shift of the metric perturbation relative to that of the aether. 
h This is not yet an observational constraint since it has not been shown that the growing 
modes are not stabilized before their effects become apparent. 
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case, if we do make this demand, then in the case of small Cj it implies 
C13 < 0, ci/cu < 1, and ci 23 /ci 4 < 1. 

On the other hand, if gravitational waves propagate subluminally rela- 
tive to the "speed of light" for matter, then matter can emit gravitational 
Cerenkov radiation. Using this phenomenon, a very tight constraint on the 
difference between the maximum speed of high energy cosmic rays and that 
of gravitational waves was derived by Moore and Nelson 27 . For cosmic rays 
of galactic origin, the constraint is Ac/c < 2 x 1CP 15 , while for extragalactic 
cosmic rays it is Ac/c < 2 x 10~ 19 . Cerenkov radiation in the additional, 
aether-metric modes has not been examined. Constraints could conceivably 
eventually be obtained using these processes. 

The requirement that all the modes propagate on the light cone of g a b 
is satisfied if and only if C4 = 0, c 3 = — ci, and c 2 = Ci/(1 — 2ci) . 



4.2. Astrophysical radiation 

Since there are three additional modes, as well as a modified speed for the 
usual gravitational waves, one expects that the energy loss rate for orbiting 
compact binaries will be affected. Not only are there additional modes to 
carry energy, but also lower multipolc moments may act as sources, for ex- 
ample there may well be a monopole moment generating the spin-0, "trace" 
mode. These phenomena have not yet been studied. Agreement with ob- 
servations of binary pulsars should yield constraints on the a parameters. 

Another potential source of constraints is the primordial perturbation 
spectrum in cosmology. Lim 14 has begun a study of this. 



5. Newtonian limit and PPN parameters 

Carroll and Lim 13 have examined the Newtonian limit. They adopt the 
ansatz of a static metric, with the aether vector parallel to the timelike 
Killing vector. Restricting to the linearized field equations they recover 
the Poisson equation V 2 J7jv — 4:TrG^p m for the gravitational potential Un, 
where p m is the usual matter energy density and 

G » - db- (10) 

where G is the parameter appearing in the action (1). Actually C4 was set 
to zero in Ref. 13, but it can be restored without calculation by using the 
fact 18 that in spherical symmetry the effect of the C4 term can be generated 
by the replacements c\ — > C14 and C3 — ► C3 — C4. 
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The PPN formalism 31 ' 32 can be applied to iE-theory since it is a metric 
theory, at least in the approximation (which is observationally known to 
be very accurate) that the matter is minimally coupled in the usual way to 
the metric. 1 In the general setting there are ten PPN parameters, but five 
of them vanish automatically in any theory that is based on an invariant 
action principle, so we need consider only the five remaining parameters, 

The two PPN parameters f3 and 7, known as the Eddington-Robertson- 
Schiff (ERS) parameters, are defined by the PPN expansion for the metric 
coefficients, 



where Un is the Newtonian gravitational potential. Thus [3 controls the 
non-linearity and 7 the space curvature due to gravity. In general relativity 
(3 = 7 = 1. The field equations are apparently too complicated to solve an- 
alytically, even assuming static, spherical symmetry. By numerical solution 
of the weak field equations in a 1/r expansion it was found 18 that the metric 
takes the form (11,12), with Un is proportional to 1/r, consistent with the 
Newtonian limit. Moreover, the ERS parameters take precisely the same 
values as in general relativity in the generic case C123 7^ 0. (For the special 
cases with C123 = see Ref. 18.) Hence the theory is indistinguishable from 
GR at the static, post-Newtonian level. 

To expose the post-Newtonian differences from GR it is necessary to 
examine the remaining PPN parameters, but these have not yet been com- 
puted for the 7E-theory The parameter £ is related to preferred location 
effects, and likely vanishes in iE-theory. The parameters are related to 
preferred frame effects and almost surely do not vanish. They will presum- 
ably arise when motion of the gravitating system relative to the asymptotic 
aether frame is allowed for. j 

One can hazard a guess based on the PPN parameters that were cal- 
culated for the vector-tensor theory without the unit constraint. 31 In that 
case, f3 and 7 are also unity in the case that corresponds most closely to 
jE-theory, namely u = in the vector-tensor parameters of Ref. 31 together 



'Strictly speaking, if the matter couples also to the aether, producing local Lorentz vio- 
lating effects in the matter sector, then the PPN formalism must be modified to allow 
for dependence of the PPN parameters on the composition of the matter sources. 
■>A recent preprint of Sudarsky and Zloshchasticv 33 addresses this point. 



.goo = 1 - 2U N + 2pU N + ■■■ 
9ij = (1 + 2~fU N H )6ij 



(11) 
(12) 
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with C4 = in our parameters 18 . If this agreement persists, one would have 
in this case £ = a\ = 0, but «2 ^ generically. However, in the special 
case when C13 = 0, for which the spin-2 waves propagate on the light cone 
of g a bi the vector-tensor parameters satisfy t = n, and the result of Rcf. 31 
yields «2 = 0. Thus perhaps in this special case all the PPN parameters of 
Einstein- iEther theory agree with those of GR. 

Current limits on a\ and a-i are of order 1CT 4 and 1CT 7 respectively, so 
constraints of this order on the parameters of iE-theory might be expected, 
at least for generic parameters. 

6. Energy 

6.1. Total energy 

As discussed in section PPN, the Newtonian potential satisfies the Poisson 
equation with source term 87rGNPm- This implies that in terms of the met- 
ric coefficient g o ~ 1 — r /r, the source mass is given by m = r /2G^. 
Assuming the source couples minimally to the metric g a b, this mass corre- 
sponds to an energy m (since the speed of light defined by g a b is one in our 
units) . Accordingly, one can infer that the energy of any isolated gravitat- 
ing system is given by the same formula. This differs from the ADM mass 
ro/2G that one would have inferred from the action (1). The reason is that 
the aether stress tensor (7) adds 1/r terms to the Newtonian field equation. 

Another path to the same conclusion is to examine the energy via its 
definition as the value of the Hamiltonian that generates asymptotic time 
translations. Using Wald's Noether charge formalism, Foster 28 showed that 
when one takes into account the falloff behavior of the fields at spatial 
infinity, the aether contribution to the energy flux integral takes the form 



where K mn a b is the tensor defined in (2). An equivalent result was found by 
Eling 29 using the Einstein pseudotensor. It was shown in Ref. 18 that for a 
spherically symmetric, static solution the line element has the asymptotic 
form 




(13) 



ds 2 = (1 - ^ + • • • ) dt 2 - (1 + ^ + • • • ) dr 2 - rW (14) 
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The r-component of the aether starts at 0(l/r 2 ), hence does not contribute 
to (13). Using this form in (13) one finds for the total energy 

Using the relation (10) between G and Gn, this can be re-expressed as 
r /2Gw Thus the total energy is related to the ADM mass Madm = r o/2G 
by a constant rescaling. 



6.2. Positivity 

For coefficients such that 2 — C14 is positive, positivity of the energy 
is thus equivalent to positivity of the ADM energy. The usual positive 
energy theorem 30 for GR assumes the dominant energy condition holds for 
the matter stress tensor, and proves that the total energy-momentum 4- 
vector of the spacetime is future timclikc. The aether stress tensor (7) does 
not appear to satisfy the dominant energy condition (for any choice of the 
Cj), so the proof does not go through as usual. Nevertheless, as discussed 
below, the energy of the linearized theory is positive for certain ranges of 
Cj. Perhaps a total divergence term that leaves the energy unchanged must 
be removed before positivity can be seen. Also, since the asymptotic value 
of the unit vector selects a preferred frame, it might be that the energy 
is always positive only in that particular frame. We can make no definite 
statement about the non-linear energy at this time, based on general formal 
grounds. 



6.2.1. Linearized wave energy 

It is useful to examine the linearized theory to begin with. The energy 
density of the various wave modes has been found 20 using the Einstein or 
Weinberg pseudotensors, averaging over oscillations to arrive at a constant 
average energy density for each mode. The energy density for the transverse 
traceless metric mode is always positive. For the vector modes it is positive 
provided (2ci — c\ + c§)/(l — C13) > 0, and for the trace mode it is positive 
provided Cu(2 — C14) > 0. For small a these energy positivity conditions 
reduce to c\ > and C14 > 0, respectively. 

In the Maxwell-like case C13 = c 2 = c 4 = the linearized energy positiv- 
ity requirement reduces to < c\ < 2. The negative energy configurations 
discussed in section 3 do not show up in the linearized limit. Their energy 
density is proportional to — (Vu ) 2 , which is quartic in the perturbation u % . 
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6.2.2. Gravitational decoupling limit 

The linearized waves are coupled metric-aether modes. A simpler limit to 
consider is a decoupling limit in which gravity is turned off. To access this 
limit formally we can let G and Cj tend to zero, while the ratios c- = Ci/G 
are held fixed. If the metric is expanded as g = n + \fGh in the action 
(1), and the limit G, C{ — > is taken, then one is left with just the action 
for linear gravitons and a decoupled aether action where all metrics are 
replaced by r\ and all covariant derivatives are replaced by ordinary partial 
derivatives. This limit was studied by Lim 14 in the case C4 = 0. Restoring 
the C4 dependence one finds perfect agreement between his results and the 
decoupling limit of the coupled linearized waves. 

As mentioned in section 3, Clayton showed that the energy can be neg- 
ative in the Maxwell-like special case. He also claimed that this remains 
true for more general choices of the coefficients Cj. However, the case where 
only ci is non-zero corresponds in the decoupling limit to a nonlinear sigma 
model (NLSM) on the unit hyperboloid which, like all NLSM's, has a stress 
tensor satisfying the dominant energy condition. 

6.3. Summary of constraints on the parameters 

So far we have discussed constraints from requirements in the linearized 
theory of positive energy, stability (no exponentially growing modes), and 
subluminal propagation (not necessarily required). Taken together, the 
constraints of positive energy and stability in the linearized regime imply 
ci > 0, C14 > 0, and C123 > (for small parameters). These are likely 
necessary for a viable theory. The requirement of no supcrluminal propa- 
gation (which we do not feel is necessary) would additionally imply C13 < 0, 
ci < C14, and C123 < C14. 

There is plenty of parameter space in which all the linearized constraints 
one might think of imposing are satisfied. It remains an important open 
question to determine whether energy positivity can be ensured beyond 
the linearized limit. But one thing is already notable, namely that this 
provides examples of a theory of a vector field which has no standard gauge 
symmetry and yet which has only stable, positive energy modes. The key 
factor making this possible is the constraint on the norm of the vector. 

7. Stars and black holes? 

In static, spherical symmetry the aether vector must be a linear combination 
of the time-translation Killing vector and the radial vector. Solutions were 
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found in Ref. 18 for which the aether field has a radial component that 
falls off as 1/r 2 . The question arises as to what happens to this vector in 
the near field region of a star. Symmetry and regularity imply that the 
radial component vanishes at the origin of spherical symmetry, so if indeed 
a regular solution exists, u a must be parallel to the Killing vector at the 
origin and at infinity but not in between. Such solutions have recently been 
studied perturbatively by Sudarsky and Zloshchastiev 33 . We have recently 
shown numerically that solutions with no radial component of the aether 
also exist 34 . 

In the case of a black hole there is no regular origin of spherical symme- 
try but the question arises as to what happens to u a on the horizon. The 
vector u a cannot exist at a bifurcation surface (where the Killing vector 
vanishes, like the 2-sphere at the origin U = V = of Kruskal coordi- 
nates). The reason is that the Killing flow acts there as a Lorentz boost in 
the tangent space of any point on the bifurcation surface, hence would act 
non-trivially on u a . Thus u a could not be invariant under the Killing flow. 
Since u a is constrained to be a unit vector it cannot vanish, hence we infer 
that it must blow up as the bifurcation surface is approached. 

One might think that the impossibility of an invariant aether at the 
bifurcation surface implies there is no regular black hole solution in this 
theory, since regularity on the future horizon is somehow connected to reg- 
ularity at the bifurcation surface. It seems this is not necessarily the case. 
A result of Racz and Wald 35 establishes, independent of any field equa- 
tions, conditions under which a stationary spacetime with a regular Killing 
horizon can be globally extended to a spacetime with a regular bifurca- 
tion surface, and conditions under which matter fields invariant under the 
Killing symmetry can also be so extended. In spherical symmetry the con- 
ditions on the metric are met for a compact Killing horizon with constant, 
non- vanishing surface gravity, so the result of Ref. 35 indicates that an ex- 
tension to a regular bifurcation surface must exist. However, one of the 
conditions on the matter (i.e. aether) field is not met, namely, it is not 
invariant under the time reflection isometry. This is because the timelike 
vector u a obviously breaks the local time reflection symmetry. Thus the 
aether vector need not be regular at the bifurcation surface (although all 
invariants must be regular and, given the Einstein equations, the aether 
stress-tensor must remain regular in the limit that the bifurcation surface 
is approached). Hence, as far as we know, there is no argument forbidding 
regular black hole solutions. 

In fact, we have expanded the field equations about a regular, static fu- 
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ture event horizon and shown 34 that locally regular solutions exist. Using a 
shooting method, we have shown that the free parameters can be chosen so 
that the solutions extend to asymptotically flat metrics at spatial infinity. 
Alternatively, it seems an attractive idea to numerically study the spher- 
ically symmetric time-dependent collapse scenario. The collapsing matter 
could be a scalar field, but more simply it could just be a spherical aether 
wave. 



8. Cosmology 

Finally we turn to the role of the aether in cosmological models. Assuming 
Robertson- Walker (RW) symmetry, u a necessarily coincides with the 4- 
velocity of the isotropic observers, so it is entirely fixed by the metric. 
The aether field equation has but one non-trivial component, which simply 
determines the Lagrange multiplier field A. Therefore the entire aether 
stress tensor is also determined by the metric. Like any matter field, when 
the aether satisfies its equation of motion, its stress tensor is automatically 
conserved. Hence, in RW symmetry, the aether stress tensor must be a 
conserved tensor constructed entirely from the spacetime geometry. One 
such tensor is the Einstein tensor itself. Another is the stress tensor of a 
perfect fluid with equation of state p = —\p, whose energy density varies 
with the scale factor a as does the spatial curvature, i.e. as 1/a 2 . The 
aether stress tensor is just a certain combination of these two conserved 
tensors 16,13 , namely 

T$ hCT = - C -ll±^l [ Gab - l^R(g ab + 2u a u b )\ . (17) 

This is written using the conventions of Refs. 18, 19 in which the field 
equations take the form (6). 

The effect of the cosmological aether is thus to renormalize the gravi- 
tational constant and to add a perfect fluid that renormalizes the spatial 
curvature contribution to the field equations. The renormalized, cosmolog- 
ical gravitational constant is given by 13 

G 

Gcosmo ~j : 7 : 7, Wo' (18) 

1 + (ci3 + 3c 2 )/2 

Carroll and Lim 13 note that, since this is not the same as Gn, the expansion 
rate of the universe differs from what would have been expected in GR 
with the same matter content. The ratio is constrained by the observed 
primordial He abundance to satisfy |G COS mo/GN — 1| < 1/8. They assume 
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the positive energy, stability, and subluminality constraints discussed above, 
which imply G cosmo < G < Gn , so the universe would have been expanding 
more slowly than in GR. In our notation, the resulting helium abundance 
constraint can be written as 15ci + 21c2 + 7c3 + 8C4 < 2, where we have 
included the c 4 dependence omitted in Ref. 13. 
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